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The 5-wave nK scattering lengths are calculated for both the isospin 1/2 and 3/2 channels in the 
lattice QCD by using the finite size formula. We perform the calculation with Nf = 2+1 gauge 
configurations generated on 32 3 x 64 lattice using the Iwasaki gauge action and nonperturbatively 
0(a) -improved Wilson action at 1/a = 2.17 GeV. The quark masses correspond to m K = 0.30 — 
0.70 GeV. For / = 1/2, to separate the contamination from excited states, we construct a 2 x 
2 matrix of the time correlation function and diagonalize it. Here, we adopt the two kinds of 
operators, su and nK. It is found that the signs of the scattering lengths are in agreement with 
experiment, namely attraction in / = 1/2 and repulsion in / = 3/2. We investigate the quark-mass 
dependence of the scattering lengths and also discuss the limitation of chiral perturbation theory. 
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1. Introduction 

The scattering length is a key quantity for understanding the basic properties of the low-energy 
interaction. Many lattice calculations of the hadron-hadron scattering lengths have been reported in 
past years. Most of them, however, do not treat the scattering system with attractive interactions due 
to the computational cost. Handling the attractive interaction would be indispensable in scattering 
studies of the future. 

Here, we focus on the 5-wave kK system. This system has two isospin channels (7 = 3/2, 1 /2). 
The low-energy interaction is repulsive (attractive) for I = 3/2 (I = 1/2). In addition, existence 
of a broad resonance is suggested in I = 1/2. Until now, three studies in lattice QCD have been 
reported [Jl], ||, ^J. The first study was performed by Miao et al. [jl|]. They calculated the 1 = 3/2 
scattering length within the quenched approximation. The first calculation with dynamical quarks 
(Nf = 2 + 1) was reported by the NPLQCD Collaboration [§]. They calculated the 1 = 3/2 scatter- 
ing length for m n = 0.3 — 0.6 GeV. They further determined the low energy constants in the SU (3) 
chiral perturbation theory (ChPT) and evaluated the 1 = 1/2 scattering length by using ChPT. The 
first direct calculation on I = 1/2 has been done by Nagata et al. [Q]. They, however, used the 
quenched approximation, and ignored the effect of ghost mesons which is nonnegligible in I = 1/2. 
In addition, their results do not reproduce the repulsive interaction for I = 3/2 at their simulation 
points. The reliability of their calculations remains controversial. In conclusion, no satisfactory 
direct calculation for I = 1/2 has been carried out. 

In the present work, we calculate the 5-wave nK scattering lengths for both the isospin chan- 
nels. We use a technique with a fixed kaon sink operator to reduce the computational cost of the 
calculation of the tcK for I = 1/2. To separate the contamination from excited states for I = 1/2, 
we construct a 2 x 2 matrix of the time correlation function and diagonalize it. After obtaining 
the scattering length at each simulation points, we investigate the quark-mass dependence and also 
discuss the limitation of &(p 4 ) SU (3) ChPT. All calculations of this work have been done on the 
super parrarell computers, PACS-CS and T2K-Tsukuba, at the University of Tsukuba. 



2. Details of simulation 

The S-wave 7tK scattering length is defined by 



an = lim tan &)(&)/& . (2.1) 

k->0 



k is the scattering momentum related to the total energy by E = ^Jrn^ + k 2 + y m 2 K + k 2 . 5q (k) is 
the S-wave scattering phase shift and can be evaluated by the Liischer's finite size formula [|J], 

where the zeta function J?oo is an analytic continuation of 

2ao(s,n) = £ * _ v (2.3) 
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defined for Re(s) > 3/2. In the case of attractive interaction, k 2 on the lowest state has a negative 
value, so k is pure imaginary. 8o(k) at the unphysical k is no longer physical scattering phase 
shift. £?oo(l,k 2 /(2n/L) 2 ), however, have a real value even for this case, so ta.n8o(k)/k obtained 
by Eq.( |2.2| ) is also real. If \k 2 \ is enough small, we can regard tanSo(k)/k as the physical scattering 
length at the nk threshold (k = 0) . 

For 1 = 3/2, one can extract E from the time correlation function 

G(0 = (0| K + ( h )% + (t) {W K+ {t + \)W n+ {t^ |0> •e" I -("- 1 ) , (2.4) 

where K + = sy$u, n + = —dy=,u, and W%+, W n + are the wall-source operators for the corresponding 
mesons. The time slice of the kaon source is shifted from that of the pion source to to avoid the 
Fierz mixing of the wall-source operators [g]. The time slice of the kaon sink operator t\ is fixed 
as t\ S> t. The exponential factor ^ n ^~ t ) is introduced to drop the unnecessary f-dependence 
appearing due to the fixed t\ . 

For / = 1/2, the existence of the K resonance are suggested in the low energy, and then it 
might be necessary to separate the contamination from the excited states. For this purpose, we use 
the two types of operators and Cl\ (Q.o and Qi), 

Qo(0 = ^(^ + (fi)7r O (0-V2^ O (^i)^ + (0)-e' M ^ ?) , 
Qi(t) = K(t) , 

Q (t ) = -j= (w K+ (t + l)W n o (to) -V2Wgo (to + l)W n+ (to)) , 

Qi(<b) = W K (ta + l), (2.5) 

where, K° = sy$d, n° = -^(Hy^u — dy^d), k = su, and W K o, W n o, W K are the wall-source operators 
for the corresponding mesons. The exponential factor in Q.o(t) is introduced by the same reason as 
for / = 3/2. We construct the 2x2 matrix of the time correlation function, 

Gij(t) = <0| Qi(t) &)(to) |0) (I, 7 = 0,1 ) , (2.6) 

and with a reference time t^ we extract the energy of the ground state by the diagonalization of 
G- l (t R ) G(t) §]. 

The calculations are carried out with Nj • = 2 + 1 full QCD configurations generated by the 
PACS-CS Collaboration ^ using the Iwasaki gauge action at j3 = 1.90 and nonperturbatively 
& (a) -improved Wilson quark action with Csw = 1-715 on 32 3 x 64 lattice. The quark propagators 
in this work are calculated with the same quark action. The corresponding lattice cutoff is l/a = 
2.176(31) GeV (a = 0.0907(13) fm) and the spatial extent of the lattice is La = 2.902(41) fm. 
The quark mass parameters and corresponding hadron masses are listed in Table [l| The Dirichlet 
(periodic) boundary condition are imposed to the temporal direction (spatial directions) in the quark 
propagators. The coulomb gauge fixing is employed for the use of the wall source. The time slice 
of the source is to = 12 (to + 1 = 13) for the 71 operator (K and K operators) and the fixed sink 
slice is t\ = 53 for K operator. We adopt tR = 18 as the reference time for the diagonalization for 
/= 1/2. The statistical errors are evaluated by the jackknife analysis with a binsize of 250 MD 
time. Here, the MD time is the number of trajectories multiplied by the trajectory length T, and 
X = 0.25 (T = 0.5) for K ud = 0.13770 (others). 
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Kud K s m n [GeV] m K [GeV] Af conf 

0.13770 0.13640 0.2971(29) 0.5935 (17) 800 

0.13754 0.13640 0.4101(27) 0.6362 (19) 450 

0.13727 0.13640 0.5707(16) 0.7135 (14) 400 

0.13700 0.13640 0.7032(11) 0.79086(97) 400 

Table 1: The quark mass parameters and corresponding hadron masses. 



3. Numerical results 

In Fig. [T], we show the time correlation functions for both the channels at m n ~ 0.30 GeV as an 
example. In the right panel, the absolute values of each component Gij(t) in / = 1/2 are presented. 
The open (filled) symbols represent the (off-)diagonal elements of Gij(t) and the signs are positive 
(negative). We find that the signals of G(t) are very clean even in I = 1/2. 
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Figure 1: (a) G{t) for / = 3/2 and(b) |G, ; (?)| for/ = 3/2 at m K ~ 0.30 GeV. 



To see effects of the contaminations from excited states for the 1=1/2 nK system, we consider 
the ratios of Goo(f) and EV [ G 1 G ( f ) ]o> the lowest eigenvalue of G 1 (t R ) G(t), to the free nK 
propagator, 

Goo(^) 



(3.1) 

In the left panel of Fig. |2[ Ro(t) (open symbols) and Do(t) (filled symbols) at m n ~ 0.30 GeV are 
plotted. The difference of the two ratios is small. This means that the contamination from the 
excited states is negligible at this quark mass parameter. On the other hand, the right panel of || 
shows that the contamination is not so small at m n ~ 0.70 GeV The diagonalization significantly 




-{m n +m K )-(t-t R ) 



D (t) = EV[G- l {t R )G(t) 



x -(m n +m K )-(t-t R ) 
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Figure 2: R (t) and D (t) at (a) m % ~ 0.30 GeV and (b) 0.70 GeV for / = 1 /2. 



changes the behavior of the ratio, because the 7r7T-type operator (On) has a large overlap with the 
excited states. Therefore, separating the contamination is indispensable for the heavy quark masses. 

k 2 and tan 8o(k)/k on the lowest state is shown in Table. ^[ For 7 = 3/2(7=1 /2), A: 2 is positive 
(negative), so we confirm the interaction is repulsive (attractive). The scattering length is defined 
as the constant term in the ^-expansion of tan 5o (&)/&. Therefore, tan<5o(&)A can be regarded as 
the scattering length if \k 2 \ is small enough to neglect the ff(k 2 ) term in the expansion. \k 2 \ for 
7 = 1/2 is, however, not so small in the heavy quark mass region. We especially find an extreme 
situation in m % > 0.41 GeV. Due to the strong attraction, tan 8q (k)/k changes the sign and we get 
tan5o(it) ~ — i. This fact suggests the appearance of an unphysical bound state of the kK system in 
m n > 0.41 GeV. We cannot use tan So(k)/k near the bound state for the extrapolation toward the kK 
threshold because the analytical structure of tan5o(&)/& is not clear. In the following discussion, 
we assume that \k 2 \ is small enough that tan $)(&)/& reflect information at the %K threshold for 
all m n of 7 = 3/2 and m K < 0.41 GeV of 7 = 1 /2, and adopt them as the scattering lengths. The 
validity of this assumption must be investigated by studying the L-dependence of tan5o(&)/& in the 
future. 





7 


= 3/2 


7 = 


1/2 


m n [GeV] 


k 2 [GeV 2 ] 


tan<5 (£)A [fm] 


k 2 [GeV 2 ] 


tan $>(*)/* [fm] 


0.30 


0.00326(30) 


-0.141(11) 


-0.00596(71) 


0.431(76) 


0.41 


0.00404(82) 


-0.170(29) 


-0.0100 (25) 


1.11 (75) 


0.57 


0.0029 (10) 


-0.130(38) 


-0.0162 (32) 


-13. (120) 


0.70 


0.00328(39) 


-0.142(15) 


-0.060 (12) 


- 0.87 (11) 



Table 2: k 2 and tan b\)(k)/k on the lowest state for 1 = 3/2 and / = 1/2. 
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We extrapolate the scattering lengths toward the physical point. For this purpose, we employ 
the formula predicted by G(p 4 ) St/ (3) ChPT ||. To improve the ChPT fit, we also include the 
data of TtTt(I = 2), which is same scattering system as TtK(I = 3/2) except for the replacement of 
the light and strange quarks. In the &(p 4 ) SU (3) ChPT, a can be described as 



a. 



a. 



(7t7tJ=2) 
Jq 



(kK, 1=3/2) 


(tcK, 1=1/2) 




\6nF 2 

SnF 2 

SnF 2 



16 

~ l + T 2 

, 16 

l + Y 2 



m 



l-LQl) + ± (ml + 2m 2 K ) ■ + % { ™^ %) 



1 



m n m K -L(fi) + ^im 2 n + 2m\) -L 4 (fi) + x^'-V 2 ) 



2 + 



16 

F 2 



m K m K -L(ju) - (ml + 2m\) ■ L 4 (n) +x ^=l/2) 



(3-2) 



where \i n ^ is the reduced mass of TtK, F is the decay constant in the chiral limit, and L4, L = 
2L\ + 2L2 +L3 — 2L4 — L5/2 + 2L(, + Lg are the low energy constants defined in Ref . [^] at scale jj, . 
^(xn:, 1=2), (n:K, 1=3 /2),(xK, 1=1/2) ^ w& ^jjQyyn functions with chiral logarithm terms and the 
explicit forms can be found in the above references. 

The fitting results of a^ n,I ~ 2 ^ jm % , a^' 1-3 ^ 2 ^ / ' \i%k and a'^ K,1 ~ X ^ T) j \l n K are plotted as a func- 
tion of m 2 K in Fig. ||[ The filled symbols represent the data used in the fit. The fit with the data in 
m K > 0.57 GeV significantly increases % /N&t, so we only use the data in m % < 0.41 GeV. The dot- 
ted lines are the fitting lines. The fit parameters and m n aQ at the physical point (m n = 0.140 GeV, 
mK = 0.494 GeV) are also listed in Table Bl where the renormalization scale is set to jl = 0.770 
GeV. The numerical results in m % < 0.41 GeV are described by the ChPT within the statistical er- 
rors. 10 3 -L 4 is consistent with -0.04(10) evaluated by PACS-CS Collaboration 0. The scattering 
lengths for kk(I = 2) and TtK (I = 3/2) are consistent with the previous lattice studies [|, 10]. Our 
measured scattering length for TtK(I = 1/2) is also consistent with the value evaluated from the 
data of TtK(I = 3/2) by using ChPT [0]. 

Our evaluation on the scattering length, however, has some issues to be solved. First, the 
statistical errors are not so small. The simple solution is to add calculations with different source 
points and improve the statistics. Second, the behavior near the chiral limit is strongly affected by 
the chiral logarithm term, so giving an evaluation without the long chiral extrapolation is desirable. 
For this purpose, the calculations in m n < 0.30 GeV are now in progress. Third, tan()o(&)/& in 
the low-momentum limit must be evaluated by systematic studies with the different volumes and 
boundary conditions. These are important issues in the future. 



X 2 /N di 


F [GeV] 10 3 -L 


10 3 -L4 m^a^ 71,1 2 ^ m w OQ K,L 3 ^ 


(nK,I=\/2) 

m„a y 


0.3 


0.114(26) 0.6(18) 


-1.4(19) -0.037(26) -0.051(38) 


0.14(10) 



Table 3: The fit parameters and m^ao at the physical point (m K = 0. 140 GeV, mjt = 0.494 GeV) 



4. Conclusion 

Direct lattice QCD computation of the S-wave scattering length of TtK (7=1/2 and / = 3/2) 
systems have been performed. The results have reproduced the correct signs of the scattering 
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(a) Repulsive channels 
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Figure 3: m^.-dependence of a /(m n ,/j. nK ) and fit curves by the &{p A ) St/ (3) ChPT. 



lengths for the first time and therefore have confirmed that the interaction is attractive (repulsive) 
in I = 1/2 (7 = 3/2). We have found that the attraction in the nK(I = 1 /2) system becomes 
stronger at m n > 0.41 GeV, and there the sign of the scattering length becomes negative. We have 
compared the m^-dependencies of the scattering lengths with those predicted by the 0(p 4 ) SU (3) 
ChPT. The data of nK(I=\ /2 and / = 3 /2), = 2) for m n < 0.41 GeV have been used in the 
fit. We have confirmed that the numerical results in m % < 0.41 GeV are described by the ChPT. 
Our evaluation on the scattering length, however, has some issues to be solved, and those remain 
as future tasks. 
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